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Beschreibung 

SIGNALTRENNUNGSVERFAHREN UND -ANORDNUNG FOR NICHTLINEARE MISCHUNGEN UNBEKAN- 
NTER SIGNALE 

Die Erfindung betrifft ein Verfahren sowie einer Anordnung zur 
Ermittlung von Parametern eines technischen Systems 

Bei einer Vielkanal-Obertragung und einem Vielkanal-Empf ang von 
Signalen tritt haufig eine Interferenz beispielsweise zwischen 
den Signalen/Bildern auf. Ein typisches Beispiel ist dabei eine 
Mischung eines Sprachsignals mit Rauschen, was bei der 
Telekommunikation und bei Videokonf erenzen ein groftes Problem 
darstellen kann. Die vorliegende Erfindung betrifft daher das 
Gebiet der Signaltrennung, urn beispielsweise ein ursprungliches 
Sprachsignal wiederzugewinnen . 

Typische bekannte Techniken zur Trennung der Quellensignale 
ausgehend von Mischsignalen basieren auf einer zeitlichen 
Mittelung oder einer Filterung der Signale. Dies hat indessen 
Nachteile bezuglich des Rechenauf wands zur Folge. 

Es sind auch Verfahren auf der Grundlage der sogenannten Blind 
Channel Equalization (Signalentzerrung ohne Vorkenntnisse uber 
den Ubertragungskanal) bekannt, aber diese Verfahren benotigen 
immer eine gewisse Kenntnis uber die Quellsignale, wie 
beispielsweise eine Kenntnis uber ihre statistische Verteilung. 

Das Problem der Signaltrennung tritt beispielsweise auch auf, 
wenn zwei Sprecher in zwei entfernt von ihnen stehende 
Mikrofone sprechen, so dafi jedes Mikrofon ein Gemisch der von 
den zwei Sprechern gesprochenen Signale auf nimmt . Somit besteht 
das Problem, das Signalgemisch, also eine Menge uberlagerter 
Eingangssignale wieder zu trennen. Aus L.Molgedey, H.G. 
Schuster, Reparation of a Mixture of Independent Signals using 
Time-Delayed Correlations", Phys. Ref. Lett. 72, 3634 (1994) 
ist dabei das folgende Verfahren bekannt: Das Problem, n 
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uberiagerte und korrelierte Quellensignale (Eingangssignale) zu 
trennen und gleichzeitige die Mischungskoeff izienten der 
Quellenstarken zu bestimmen, lafit sich auf ein Eigenwertproblem 
reduzieren, bei dem simultan zwei symmetrische n x n Matrizen 
diagonalisiert werden miissen. Die Matrixelemente sind meftbare 
zeitverzogerte Korrelationsf unktionen. Die Losung des 
Eigenwertproblems kann durch ein neuronales Netz erfolgen, 
wobei die Lernregel fur die lateralen inhibitorischen 
Wechselwirkungen zwischen den Neuronen durch eine Liapunov- 
Funktion bestimmt werden kann, deren Minima die (entarteten) 
Losungen des Problems liefern. 

Dieses Verfahren wurde auch bereits auf akustische 
Eingangssignale angewandt (siehe F. Ehlers, H. G. Schuster, 
„Blind Separation of convolutive mixtures and an application in 
automatic speech recognition" IEEE Trans. Signal Proc. (1997). 

Aus der DE 195 31 388 CI ist ein Signaltrennungsverf ahren und 
eine Signaltrennungseinrichtung fur nichtlineare Mischungen 
unbekannter Signale (Blind Channel) bekannt, das schematisch in 
Fig. 3 dargestellt ist. 

Dieses deutsche Patent behandelt die Separierung eines 
Signalgemisches, bestehend aus der nichtlinearen Uberlagerung 
von M unbekannten Quellsignalen XI, X2, wobei N (N > M) 
unterschiedliche Mischungen der M Quellsignale XI, X2 inklusive 
eines eventuellen Storsignals einer Signalauswerteeinrichtung 
zugefuhrt werden, die durch eine statistische Analyse der 
Signale die nichtlinearen Obertragungsf aktoren bestimmt und mit 
diesen errechneten Faktoren die Mischung riickgangig macht, so 
daB die N Ausgange der Signaltrennungseinrichtung moglichst 
naherungsweise die M Quellsignale ohne Oberlagerungen 
enthalten. Dadurch wird eine Behandlung nichtlinearer Gemische 
moglich, wobei nichtlinear bedeutet, dafi die Quellsignale XI, 
X2 durch ein unbekanntes nichtlineares System G gemischt 
werden. Das unbekannte System G wird durch eine sogenannte 
Volterra-Reihe beschrieben, und die Signaltrennungseinrichtung 
G-l bestimmt die Koeff izienten der Volterra-Reihe, Mit deren 
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Kenntnis ist eine Entmischung des Signalgemischs moglich. 
Aufierdem konnen die Koef f izienten zu weiteren Analysen zur 
Orts- oder Geschwindigkeitsbestimmung der Signalquellen benutzt 
werden. 

Das aus dieser Druckschrift bekannte Verfahren besteht dabei im 
wesentlichen aus zwei Schritten: 

• Erstens werden die durch den wahlbaren Grad der 
Nichtlinearitat bei der Mischung eindeutig bestimmten 
nichtlinearen Gleichungen flir ein gleitendes Zeitfenster 
gelost und die Losungen werden iiber die Zeit gemittelt. Diese 
zeitliche Mittelung stellt einen Hauptnachteil dieser 
bekannten Technik dar, da sie den Rechenauf wand und 
gleichzeitig die Zeitdauer fur die Berechnung erhoht. 

• Zweitens wird aus einer genugend groften Anzahl von 
unterschiedlichen Kumulanten der geschatzten Ausgangssignale 
gebildetes Potential minimiert , wobei die zur Berechnung des 
Potentials notigen Werte aus einem gleitenden Zeitfenster 
einer wahlbaren Lange stammen. Dabei ist vorausgesetzt , daft 
sich das Mischungssystem so langsam andert, daft diese 
Anderung bei der Berechnung der gesuchten Mischungsf aktoren 
vernachlassigt werden kann. Gemaft diesem deutschen Patent 
wird zur Ausfuhrung des zweiten genannten Schrittes eine 
Kostenf unktion konstruiert, die minimiert wird. Wenn das 
globale Minimum erreicht ist, hat man die optimalen Werte, in 
diesem Fall der Ubertragungsf aktoren, gefunden. 

Hinsichtlich des zeitlichen Aufwands sowie des Rechenauf wands 
ist das in der DE 195 31 388 CI beschriebene Verfahren aufgrund 
der Ausfuhrung einer zeitlichen Mittelung am Abschluft des 
ersten, oben genannten Verf ahrensschrittes nachteilig. 

Die vorliegende Erfindung hat daher zur Aufgabe, ein Verfahren 
und eine Anordnung bereitzustellen, die die Trennung 
uberlagerter., statistisch voneinander unabhangiger akustischer 
Signale mit verringertem Rechenauf wand ermoglicht. 
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Diese Aufgabe wird durch das Verfahren sowie durch die 
Anordnung mit den Merkmalen gemali den unabhangigen Anspruchen 
gelost . 

Bei einem Verfahren zur Ermittlung von Parametern eines 
technischen Systems, mit dem Ausgangssignale aus einer Menge 
iiberlagerter , statistisch voneinander unabhangiger 
Eingangssignale ermittelt werden kbnnen, werden die Parameter 
derart ermittelt, dafi die statistische Unabhangigkeit der 
Ausgangssignale maximiert wird. 

Eine Anordnung zur Ermittlung von Parametern eines technischen 
Systems, mit dem Ausgangssignale aus einer Menge iiberlagerter , 
statistisch voneinander unabhangiger Eingangssignale ermittelt 
werden konnen, weist einen Prozessor auf, der derart 
eingerichtet ist, daft die Parameter derart ermittelt werden 
konnen, dali die statistische Unabhangigkeit der Ausgangssignale 
maximiert wird. 

Vorteilhafte Weiterbildungen der Erfindung ergeben sich aus den 
abhangigen Anspruchen. 

Die Parameter werden bevorzugt in einem iterativen Verfahren 
ermittelt. 

In einer weiteren Ausgestaltung sind die Parameter Elemente 
einer Entmischmatrix, mit der die Menge der uberlagerten 
Eingangssignale multipliziert oder auch gefaltet wird, wodurch 
die Ausgangssignale gebildet werden. 

Die Optimierung der Parameter der Entmischmatrix wird bevorzugt 
durch die folgenden Schritte erhalten: 

- Wiederholung einer zeitverzogerten Dekorrelationsberechnung 
zur Ermittlung der Eigenwerte der Entmischmatrix, 
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- Ermittlung der Eigenwerte der Entmischmatrix, fiir die 
Kreuzkorrelationen einen minimalen Wert annehmen, und 

- Ausfiihrung einer Kumulantenminimierung, wobei als Startwerte 
fur die Kumulantenminimierung die im vorherigen Schritt 
ermittelten Eigenwerte verwendet werden. 

Die Kumulantenminimierung kann beispielsweise durch Trainieren 
eines neuronalen Netzes, aber auch jegliche andere bekannte 
Minimierungstechnik, wie beispielsweise Gradientenabstieg oder 
Monte-Carlo-Simulationen verwendet werden. 

In einer Weiterbildung wird bei der Optimierung der Parameter 
der Entmischmatrix wenigstens ein Diagonalparameter der 
Entmischmatrix auf einen vorgegebenen Wert gesetzt, womit , um 
eine Stabilitat des Minimierungsvorgangs hin zu einem globalen 
Minimum zu gewahrleisten . 

Die Entmischmatrix wird bevorzugt auf eine finite Impulsantwort 
begrenzt, d.h. es wird ein FIR-Filter (Finite _Impulse Response) 
eingesetzt zur Bildung der einzelnen Komponenten der 
Entmischmatrix. Der FIR-Filter kann sowohl ein kausales FIR- 
Filter oder auch ein nicht-kausales FIR-Filter sein. 

Ferner wird die Entmischmatrix wahrend der 

Kumulantenminimierung bevorzugt durch Projektion in einen 
Einheitskreis stabilisiert . 

Die Weiterbildungen gelten sowohl fur das Verfahren als auch 
fiir die Anordnung, wobei bei der Anordnung jeweils der 
Prozessor derart eingerichtet ist, dafi der entsprechende 
Verf ahrensschritt durchfiihrbar ist oder durchgefiihrt wird. 

Die Erfindung sowie deren Weiterbildungen konnen vorteilhaft 
eingesetzt werden zur Trennung uberlagerter , statistisch 
voneinander unabhangiger Eingangssignale, insbesondere 
akustischer Eingangssignale. 
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Das Verfahren sowie die Anordnung sind fur eine beliebige 
Anzahl von Eingangssignalen anwendbar. 

Weitere Vorteile, Merkmale und Eigenschaf ten der vorliegenden 
Erfindung werden nunmehr bezugnehmend auf die beiliegenden 
Figuren der Zeichnungen anhand eines Ausf uhrungsbeispiels naher 
erlautert . 

Fig. 1 zeigt die Anwendung eines Systems zur Trennung 
iiberlagerter , statistisch voneinander unabhangiger 
akustischer Signale gemaft dem Ausf uhrungsbeispiel , 

Fig. 2 zeigt eine symbolische Darstellung des Systems von 
Fig. 1, und 

Fig. 3 zeigt eine aus dem Stand der Technik (DE 195 31 388 
CI) bekannte Signaltrennungseinrichtung fur nichtlineare 
Mischungen unbekannter Signale. 

Zur Wiedergewinnung des ursprunglichen Sprachsignals wird aus 
einer Mischung von Signalen die statistische Unabhangigkeit 
zwischen den Quellensignalen (ursprungliches Sprachsignal und 
Rauschen) , im weiteren auch als Eingangssignale bezeichnet, 
ausgenutzt und der inverse Vorgang des dynamischen Systems, der 
die Mischung der Signale ergeben hat, wird im wesentlichen 
naherungsweise trainiert (gelernt) . Zwei verschiedene 
Mischungen des Sprachsignals bzw. des Rauschsignals werden 
beispielsweise durch zwei Mikrofone 1, 2 (vgl. Fig.l) erhalten, 
die voneinander beabstandet sind und/oder in entgegengesetzten 
Richtungen ausgerichtet sind. Bei dem Verfahren wird die 
sogenannte zeitverzogerte Dekorrelationstechnik (TDD, time 
delayed decorrelation) verwendet, urn die Lernphase zu 
initiieren, d.h. urn Startwerte fur die Lernphase zu ermitteln 
und vorzugeben, wodurch der Berechnungsauf wand fur eine im 
weiteren beschriebenen Kumulantenminimierung verringert werden 
kann und die Gefahr lokaler Minima verringert werden kann. 
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Fig.l zeigt zwei Mikrofone 1, 2, die ein erstes Eingangssignal 
Zl(t) und ein zweites Eingangssignal Z2(t) aufnehmen. Diese 
Eingangssignale Zl(t) und Z2(t) konnen untereinander wiederum 
jeweils mit Rauschen vermischt sein, was durch eine Mischmatrix 
S (siehe Bezugszeichen 3) symbolisch in Fig.l dargestellt ist. 
Nach dem Empfang bzw. der Obertragung wird eine Menge XI (t) und 
X2(t) uberlagerter, statistisch voneinander unabhangiger 
Eingangssignale Zl(t) und Z2(t) erhalten. Diese Signale werden 
in eine Berechnungseinheit 4 eingegeben, in der im wesentlichen 
zwei Schritte ausgefiihrt werden, die symbolisch durch eine 
Berechnungseinheit B (Bezugszeichen 6) fur den ersten Schritt 
sowie ein neuronales Netzwerk 5 fur den zweiten Schritt 
dargestellt ist - 

Durch die Berechnungseinheit 4 werden zwei Ausgangssignale 
Yl(t) bzw. Y2(t) ermittelt, die bei optimaler Einstelllung der 
Parameter in der Berechnungseinheit 4 naherungsweise gleich den 
Eingangssignalen Zl(t) bzw. Z2(t) sind. Mit anderen Worten, bei 
optimaler Einstellung der Parameter der verwendeten Matrizen in 
der Berechnungseinheit 4 erfolgt durch diese Berechnungseinheit 
4 im wesentlichen der inverse Vorgang zu dem dynamischen 
Mischvorgang, der durch die Matrix S (Bezugszeichen 3) 
symbolisch dargestellt ist. Das Ausf uhrungsbeispiel beschaftigt 
sich mit dem Optimierungsvorgang der Einstellung der Parameter 
der Entmischmatrix. 

Es werden die Parameter der Matrizen in der Berechnungseinheit 
4 derart optimiert, dafi die statistische Unabhangigkeit 
zwischen den durch den Matrizenvorgang in der 
Berechnungseinheit 4 gewonnenen Ausgangssignale Yl(t), Y2(t) 
maximiert wird. Zu diesem Zweck werden die Ausgangssignale 
Yl(t) bzw. Y2(t) zu der Berechnungseinheit 4 zuruckgef iihrt (s. 
Ruckfuhrschleif en 7 bzw. 8) . Durch ein iteratives Verfahren 
wird ermittelt, ob sich die statistische Unabhangigkeit der 
Ausgangssignale Yl(t) bzw. Y2(t) im Vergleich zu dem vorherigen 
Schritt der Iteration erhoht hat (und somit die Iteration die 
„richtige" Richtung in Richtung des globalen Minimums einer im 
weiteren beschriebenen Kostenf unktion einnimmt) oder nicht . 
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Fig. 2 zeigt eine mathematische Darstellung des Schemas von 
Fig.l, wobei der Mischvorgang 3 durch eine Matrix S(q) 
mathematisch beschrieben werden kann und der Entmischvorgang, 
der durch die Berechnungseinheit 4 erfolgen soli, durch eine 
Entmischmatrix M(q) symbolisiert wird. 

In Fig. 2 ist somit das Problem der Trennung einer sogenannten 
Multi-Channel Blind Source (Vielf ach-Kanal-Quelle ohne a- 
priori-Kenntnis) in zwei Dimensionen dargestellt. Dabei ist 
angenommen, daft das Mischsystem S(q), wobei q fur eine 
Einheitsverzogerung steht, stabil ist und gleichzeitig auch 
eine stabile Invertierung aufweist, d.h. dafi es ein 
Minimalphasensystem ist. Dariiber hinaus ist angenommen, daft die 
Eingangssignale Zl(t) und Z2(t) (beispielsweise ein Sprach- 
bzw. ein Rauschsignal) statistisch voneinander unabhangig sind 
und nicht gaussformig verteilt sind. Die Menge XI (t) und X2(t) 
der uberlagerten Eingangssignale Zl(t) und Z2(t) sind 
Eingangssignale in ein Entmischsystem mit einer Entmischmatrix 
M(q), deren Parameter (Matrixelemente) auf eine Maximierung der 
statistischen Unabhangigkeit zwischen den Ausgangssignalen 
Yl(t) und Y2(t) trainiert werden. Unter „Trainieren" ist dabei 
der gut bekannte Lernvorgang beispielsweise eines neuronalen 
Netzes bekannt, das als ein Beispiel fur eine Technik genannt 
sein soil, die statistische Unabhangigkeit zu maximieren. Dies 
erfolgt durch Minimierung einer im weiteren beschriebenen 
Kostenf unktion J(M). 

Es wird eine Kumulanten-Kostenf unktion gebildet, die die 
Diagonalkumulantenelementen der Kumulantenordnung 2-4 
minimiert : 

4 . 2 

Dcum s J(M) = ]T ^ C (l) nondiag 
i = l nondiag 



Folgende Gesichtspunkte der dynamischen Mischung durch die 
Mischmatrix S(q) sind dabei zu berucksichtigen : 



WO 00/25489 




PCT/DE99/03304 



9 

Stabilitat des Entmischsystems : 
Dies wird erreicht, wenn M(q) auf eine finite 
Impulsantantwort (FIR-Filter) beschrankt wird. Die 
Stabilitat des FIR-Systems M(q) kann daruber hinaus auch 
dadurch erhalten werden, dafi wahrend der Lernphase eine 
Projektion in den Einheitskreis erfolgt Eine 
moglicherweise vorliegende Nichtkausalitat der 
Invertierung von S(q) kann durch eine geeignete 
Zeitverschiebung (Verzogerung) des Eingangssignals X(t) 
kompensiert werden . 

Eindeutigkeit der getrennten Signale Y(t): 
Fur den Fall statischer Mischungen werden die 
ursprunglichen Quellsignale durch eine Skalierung 
wiedergewonnen. Fur den Fall einer dynamischen Entmischung 
ist die Gefahr einer Nichteindeutigkeit der getrennten 
Signale Y(t) sogar noch grofier. Es ist of f ensichtlich, daii 
fur den Fall, daft Yl(t) und Y2(t) statistisch voneinander 
unabhangig sind, auch jegliche linear gefilterte 
Modifikation dieser Signale immer noch statistisch 
unabhangig sind. Daher ist eine zusatzliche Information 
notwendig, urn die inharente Nichteindeutigkeit des 
Problems zu verringern. 

- Gaussverf ormung der Daten: 

Algorithmen auf Kumulantenbasis fur eine statische Blind 
Source-Trennung minimieren bzw. eliminieren in effektiver 
Weise Diagonalkumulanten hoherer Ordnung entsprechend den 
Ausgangssignalen Y(t). Andererseits fuhrt eine lineare 
Filterung zu einer Gaussverteilungs-Verf ormung der Daten, 
bei denen die Kumulanten hoherer Ordnung in Richtung 0 
gehen. Dies kann daher zu Randlosungen fuhren, in denen 
die Kostenf unktion ein lokales Minimum erreicht, wobei die 
erwunschte eigentliche Trennung (globales Minimum) nicht 
erfolgt. Urn diesen ungewunschten Fall zu vermeiden, wird 
die Struktur der Entmisch-Transf erf unktion 
(Entmischmatrix) M(q) einigen Beschrankungen unterworfen. 



WO 00/25489 




PCT/DE99/03304 



Urn die obengenannten Probleme zu umgehen, wird ein Ansatz 
gewahlt, daft wenigstens eines (oder auch alle) der 
Diagonalelemente auf den Einheitswert gesetzt werden: 

Mll(q) - 1 

und 

M22 (q) =1. 

Diese Annahme ist exakt, wenn die Mischelemente Sll(q) und/oder 
S22(q) ebenfalls einen Einheitswert „1" aufweisen. Sonst sei 
angenommen, daft Sll (q) und/oder S22(q) stabile Invertierungen 
aufweisen, was die Skalierung der Diagonalelemente von M(q) auf 
den Einheitswert erlaubt, Dieser Ansatz verringert wesentlich 
die Nichteindeutigkeit einer Losung und vermeidet dariiber 
hinaus in effektiver Weise die Gefahr einer ubermaftigen 
Gaussverteilungs-Verf ormung der Ausgangssignale . Auch wenn auf 
den ersten Blick die Beschrankung der Diagonalelemente von 
M(q), wie oben ausgefuhrt, sehr restriktiv wirkt, ist diese 
Annahme in der praktischen Anwendung in der Regel erfullt. Ein 
typisches Beispiel ist die Entrauschung von Sprachsignalen auf 
der Grundlage einer Aufzeichnung mit zwei Mikrofonen, wobei die 
Mikrofone voneinander raumlich getrennt sind oder ein Mikrofon 
in Richtung des Sprechers gerichtet ist, wahrend das andere 
Mikrofon in der umgekehrten Richtung gerichtet ist, so daft das 
zweite, von dem Sprecher abgewandte Signal im wesentlichen nur 
ein Rauschsignal aufnimmt. 

Der Kumulantenansatz basiert auf einer direkten Ermittlung der 
Diagonalkumulanten, wie sie in dem eingangs genannten Artikel 
von F. Ehlers und H. Schuster, „Blind Separation of convolutive 
mixtures and an application in automatic speech recognition" 
IEEE Trans. Signal Proc. (1997), ausgefuhrt ist, was aber 
inharent den Nachteil aufweist, daft die numerische Losung sehr 
aufwendig ist. Daher wird eine geeignete Initialisierung dieses 
Minimierungsverf ahrens verwendet. Urn Startwerte fur das 
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Minimierungsverf ahren zu ermitteln, wird gemaft der vorliegenden 
Erfindung die Technik der zeitverzogerten Dekorrelation (TDD) 
zur gleichzeitigen Dekorrelierung zweier verschiedener 
Zeitverzogerungen angewendet, wobei diese TDD-Technik sich auf 
ein geeignetes Matrix-Eigenwertproblem zurlickfuhren lafit- Wie 
bereits gesagt, wird diese TDD-Technik gemaft der vorliegenden 
Erfindung zur Initiierung des Diagonal (Kreuzkorrelations) - 
Kumulantenminimierungsproblems verwendet . 

Zusammenf assend laftt sich das Verfahren in die folgenden beiden 
Schritte unterteilen: 

1. Wiederholung des TDD-Verf ahrens auf Grundlage des 
Eigenwertproblems im Frequenzbereich fur verschiedene 
Verzogerungspaare und Ermittlung der Losung, bei der die 
Kreuzkorrelationsterme einen minimalen Wert aufweisen. 

2. Initiierung (Start) der Diagonal-Kumulantenminimierung auf 
Grundlage der im oben genannten Schritt ermittelten 
Startwerte ( FIR-Parameter ) . 

Im folgenden sollen noch einige Haupteigenschaf ten und Vorteile 
noch einmal zusammengef aftt werden: 

- Es ist kein a-priori-Wissen der Signaleigenschaf ten 
notwendig, mit der Ausnahme, daft eine statistische 
Unabhangigkeit gefordert wird. 

- Die Stabilitat des dynamischen Entmischsystems wird durch 
die Modulierung seiner Komponenten als FIR-Filter 
gewahrleistet . 

- Eine ubermaftige Gaussverteilung-Verf ormung wird durch den 
Ansatz vermieden, daft wenigstens eines der Elemente der 
Mischtransf er-Funktionsmatrix (Entmischmatrix) auf den 
Einheitswert gesetzt wird bzw. auf den Einheitswert 
skaliert werden kann, und 
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- da der Kumulanten-Minimierungsschritt (Schritt 2) einen 
hohen Rechenauf wand erfordert, wird der Lernalgorithmus 
beispielsweise eines neuronalen Netzes durch das TDD- 
Ver f ahren initialisiert • 



WO 00/25489 



PCT/DE99/03304 



13 

Im weiteren ist ein Programm in Matlab, Version 4 oder Version 
angegeben, mit dem das oben beschriebene Ausf uhrungsbeispiel 
auf einem Rechner realisiert werden kann: 

function [ cost , outl , out2 ] =cumulant_cost FIRa2 (par, input , pi , pll , 
p2,p22, a3, a4) ; 

% [cost, outl, out2] =cumulant_costFIRa2 (par, input, pi, pll, p2, p22, 
a3, a4 ) ; 

% cumulant cost 

% FIR representation used 

% filter function used in both directions (non_causal) 

[np,mp] =size (par ) ; 
firl=par (l:pl) ; 
firll=par (l+pl:pl+pll) ; 
fir2=par(l+pl+pll:pl+pll+p2) ; 
fir22=par (pl+pll+p2+l :mp) ; 

den=l; %FIR only 

outl- [input ( : ; 1) -filter (firl, den, input ( : ,2) ) -flipud (filter 
([0 firll ], [den] , flipud (input (:, 2) )))] ; 

%/std (input (1, : ) ) ; %dlsim 

%f ilter 

out2= [input ( : , 2) -f ilter ( fir2 , den, input ( : , 1) ) - flipud (filter 
([0 fir22 ] , [den] , flipud (input ( : , 1) ) ) ) ] ; %/std (input ( : ,2) ) ; 

%dlsim %filter 

out = [ out 1 out 2 ] ; 

%outl=outl/std (outl) ; % this scaling was not needed in examples 
in SIP98 paper 
%out2=out2/std (out2) ; 

Ld=0; % number of delays in calculating the cross-correlation 

cost3=0; 

cost4=0 ; 
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COStALLl=[] ; 
costALL2=[] ; 
ol2=outl. *out2; 
cost2=mean (ol2) A 2; 

oll2=outl . *outl • *out2 ; 
ol22=outl. *out2.*out2; 

if a3 -=1 

cost3=[mean(ol22) ] A 2 + [mean (ol22 ) ] "2; 

end 

if a4 =-1 

cost4= [mean (oll2 . *outl) -3*mean (outl . A 2) *mean (ol2) ] A 2+. . . 
[mean( (outl. A 2) . *(out2. A 2) ) -2*mean <ol2) A 2- 
mean (outl. A 2) *mean (out2 . A 2 ) ] A 2+. . . 
[mean(ol22.*out2)-3*mean(out2. A 2) *mean(ol2) ] A 2; 



end 



%cum4a= [cum4x (outl, outl, outl, outl) ] A 2; 
%cum4b= [cum4x (out2, out2, out2, out2) ] A 2; 



cost=cost2+a3*cost3+a4*cost4 ; 



cum4a-cum4b; 
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Paten tanspriiche 

1. Verfahren zur Ermittlung von Parametern eines technischen 
Systems, mit dem Ausgangssignale aus einer Menge uberlagerter, 
statistisch voneinander unabhangiger Eingangssignale ermittelt 
werden konnen, bei dem die Parameter derart ermittelt werden, 
daft die statistische Unabhangigkeit der Ausgangssignale 
maximiert wird. 

2. Verfahren nach Anspruch 1, 

bei dem die Parameter in einem iterativen Verfahren ermittelt 
werden . 

3. Verfahren nach Anspruch 1 oder 2, 

bei dem die Parameter Elemente einer Entmischmatrix sind, mit 
der die Menge der uberlagerten Eingangssignale multipliziert 
wird, wodurch die Ausgangssignale gebildet werden. 

4. Verfahren nach Anspruch 3, 

bei dem die Optimierung der Parameter der Entmischmatrix durch 
die folgenden Schritte erhalten wird: 

- Wiederholung einer zeitverzogerten Dekorrelationsberechnung 
(6) zur Ermittlung der Eigenwerte der Entmischmatrix, 

- Ermittlung der Eigenwerte der Entmischmatrix, fur die 
Kreuzkorrelationen einen minimalen Wert annehmen, und 

- Ausfiihrung einer Kumulantenminimierung (5), wobei als 
Startwerte fur die Kumulantenminimierung die im vorherigen 
Schritt ermittelten Eigenwerte verwendet werden. 

5. Verfahren nach Anspruch 4, 

bei dem die Kumulantenminimierung durch Trainieren eines 
neuronalen Netzes (5) erfolgt. 

6. Verfahren nach einem der Anspriiche 3 bis 5, 
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bei dem bei der Optimierung der Parameter der Entmischmatrix 
wenigstens ein Diagonalparameter der Entmischmatrix auf einen 
vorgegebenen Wert gesetzt wird. 

7. Verfahren nach einem der Anspriiche 3 bis 6, 

bei dem die Entmischmatrix auf eine finite Impulsantwort 
begrenzt wird. 

8. Verfahren nach einem der Anspriiche 3 bis 7, 

bei dem die Entmischmatrix wahrend der Kumulantenminimierung 
(5) durch Projektion in einen Einheitskreis stabilisiert wird. 

9. Verfahren nach einem der Anspriiche 1 bis 8, eingesetzt zur 
Trennung iiberlagerter , statistisch voneinander unabhangiger 
Eingangssignale . 

10. Verfahren nach einem der Anspriiche 1 bis 8, eingesetzt zur 
Trennung iiberlagerter , statistisch voneinander unabhangiger 
akustischer Eingangssignale. 

11. Anordnung zur Ermittlung von Parametern eines technischen 
Systems, mit dem Ausgangssignale aus einer Menge uberlagerter, 
statistisch voneinander unabhangiger Eingangssignale ermittelt 
werden konnen, mit einem Prozessor, der derart eingerichtet 
ist, daft die Parameter derart ermittelt werden konnen, daii die 
statistische Unabhangigkeit der Ausgangssignale maximiert wird, 

12. Anordnung nach Anspruch 11, 

bei der der Prozessor derart eingerichtet ist, dafi die 
Parameter in einem iterativen Verfahren ermittelt werden. 

13. Anordnung nach Anspruch 11 oder 12, 

bei der der Prozessor derart eingerichtet ist, daii die 
Parameter Elemente einer Entmischmatrix sind, mit der die Menge 
der uberlagerten Eingangssignale multipliziert wird, wodurch 
die Ausgangssignale gebildet werden. 
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14. Anordnung nach Anspruch 13, 

bei der der Prozessor derart eingerichtet ist, daft die 
Optimierung der Parameter der Entmischmatrix durch die 
folgenden Schritte erhalten wird: 

- Wiederholung einer zeitverzogerten Dekorrelationsberechnung 
(6) zur Ermittlung der Eigenwerte der Entmischmatrix, 

- Ermittlung der Eigenwerte der Entmischmatrix, fur die 
Kreuzkorrelationen einen minimalen Wert annehmen, und 

- Ausfiihrung einer Kumulantenminimierung (5), wobei als 
Startwerte fur die Kumulantenminimierung die im vorherigen 
Schritt ermittelten Eigenwerte verwendet werden. 

15. Anordnung nach Anspruch 14, 

bei der der Prozessor derart eingerichtet ist, daft die 
Kumulantenminimierung durch Trainieren eines neuronalen Netzes 
(5) erfolgt. 

16. Anordnung nach einem der Anspruche 13 bis 15, 

bei der der Prozessor derart eingerichtet ist, daft bei der 
Optimierung der Parameter der Entmischmatrix wenigstens ein 
Diagonalparameter der Entmischmatrix auf einen vorgegebenen 
Wert gesetzt wird. 

17. Anordnung nach einem der Anspruche 13 bis 16, 

bei der der Prozessor derart eingerichtet ist, dali die 
Entmischmatrix auf eine finite Impulsantwort begrenzt wird. 

18. Anordnung nach einem der Anspruche 13 bis 17, 

bei der der Prozessor derart eingerichtet ist, daft die 
Entmischmatrix wahrend der Kumulantenminimierung (5) durch 
Projektion in einen Einheitskreis stabilisiert wird. 



19. Anordnung nach einem der Anspriiche 11 bis 18, eingesetzt 
zur Trennung- uberlagerter , statistisch voneinander unabhangiger 
Eingangssignale . 
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20. Anordnung nach einem der Anspruche 11 bis 18, eingesetzt 
zur Trennung uberlagerter , statistisch voneinander unabhangige 
akustischer Eingangssignale . 
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